From magnetism to one-dimensional spin liquid in the anisotropic triangular lattice 
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We investigate the anisotropic triangular lattice that interpolates from decoupled one-dimensional 
chains to the isotropic triangular lattice and has been suggested to be relevant for various quasi- 
two-dimensional materials, such as CS2CUCI4 or k-(ET)2Cu2(CN)3, an organic material that shows 
intriguing magnetic properties. We obtain an excellent accuracy by means of a novel representation 
for the resonating valence bond wave function with both singlet and triplet pairing. This approach 
allows us to establish that the magnetic order is rapidly destroyed away from the pure triangular 
lattice and incommensurate spin correlations are short range. A non-magnetic spin liquid naturally 
emerges in a wide range of the phase diagram, with strong one-dimensional character. The relevance 
of the triplet pairing for k-(ET)2Cu2(CN)3 is also discussed. 

PACS numbers: 75.10.-b, 71.10.Pm,75.40.Mg 



When cooling down the temperature, the majority of 
materials undergo phase transitions to ordered phases, 
that break some symmetry. Examples are ubiquitous 
in nature, e.g., magnets or superconductors, and define 
the paradigm in solid state physics. In the last years, a 
great effort has been done to determine and character- 
ize new states of matter, which escape this conventional 
description. In this regard, one of the most intriguing 
case is given by the so-called spin liquids, namely in- 
sulating phases that cannot be adiabatically connected 
to any band insulators. 1] The concept of spin liquid 
was introduced by Fazekas and Anderson 2] and its 
possible connection with the low-doping regime of high- 
temperature superconductors was highlighted by Ander- 
son. ^ The standard picture of a spin liquid is given by 
the resonating valence bond (RVB) ansatz, a superposi- 
tion of configurations in which couples of spins form sin- 
glets but change partner from one configuration to the 
other. After a long period dominated by the (wrong) 
prejudice that spin liquids cannot be actually stabilized, 
today there is an increasing evidence that they can be 
obtained in both microscopic models and real materi- 
als. Spin-liquid behavior has been suggested in vari- 
ous compounds: in two-dimensional (2D) triangular lat- 
tices, Aij^l hi Kagome materials, Q and more recently in 
three-dimensional hyper-Kagome antiferromagnets. [7[ 

The 2D triangular lattice is the simplest structure in 
which the nearest-neighbor super-exchange leads to frus- 
tration. However, it is well proved that ideal Heisenberg 
spins with antiferromagnetic interactions on such a lat- 
tice display an ordered spin configuration, even for the 
spin-half case. Nevertheless, due to strong quantum 
fluctuations, the magnetic order parameter is highly re- 
duced from its classical value 0, [l3] and small perturba- 
tions may destroy long-range order and drive the system 
towards a pure spin-liquid ground state. In this sense, a 
finite on-site repulsion U (or equivalently multi-spin in- 
teractions) may stabilize a magnetically disordered phase 



close to the metal-insulator transition. 11, l2| Another 
very interesting possibility to further increase quantum 
fluctuations is to have different super-exchange couplings 
along different spatial directions. This latter case is par- 
ticularly appealing because of its connection with various 
materials, such as CS2CUCI4 and Cs2CuBr4 [l3| or a fam- 
ily of quasi- 2D organic compounds. 14 1 

In this Letter, we consider a spin-half Heisenberg 
model defined on the anisotropic triangular lattice: 
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{iJ} 



where Si — {Sf, Sf , S'f ) is the spin operator at the site i 
and indicates nearest-neighbor sites along the ai = 
(1,0) direction, while indicates nearest-neighbor 

sites along either a2 = (1/2, \/3/2) or ag = (-l/2,V5/2). 
Therefore, the model consists in one-dimensional (ID) 
chains coupled with zig-zag bonds J'. Here, we consider 
clusters with N = sites and periodic boundary condi- 
tions along Lai and La2. Recent works showed a strong 
one-dimensionalization [l5| and gapless S* = 1/2 excita- 
tions [l6[ in a wide regime of frustration J'/ J < 0.5. The 
main limitation of these results is the inaccurate descrip- 
tion of the magnetic correlations. In fact, works based 
upon series expansions 1^, 1^ showed that a magnetic 
spiral order may be present down to the ID limit, with 
almost antiparallel spins along chains. The fact that the 
ID disordered phase is unstable towards the formation of 
incommensurate magnetic order has been also suggested 
by perturbative approaches. 2^, 21 1 The situation is far 
from being clarified and more work is needed to under- 
stand the nature of the ground state. 

Before considering our variational Monte Carlo calcu- 
lations, it is useful to present exact results by the Lanczos 
method on the 6x6 cluster, see Fig. [TJ As already ob- 
tained in Ref. l22l . we find a level crossing in the ground 
state; this is due to a change, around J'/ J ~ 0.825, in the 
quantum number of the reflection symmetry. In addition 
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FIG. 1: (Color on-line) Exact energy gap for the 6x6 cluster 
as a function of the frustrating ratio J'/J. Full and empty 
circles indicate singlet states with q = (0, 0) and different 
reflection symmetry (the data show a level crossing for J'/ J ~ 
0.825). Empty and full triangles indicate triplet excitations at 
q = (7r,7r/3\/3) and q — (7r,7r/\/3), respectively. Full squares 
indicate triplet excitations at g = (37r/4, 0). 



FIG. 2: (Color on-line) Accuracy for the energy on the 6x6 
lattice, Eq and E denote the exact and the variational en- 
ergies. The WF with decoupled chains (black triangles) and 
the one with 2x1 structure and 120° order (red squares) are 
reported, see text for a detailed description. The example of 
the 2x1 structure for the sign of the nearest-neighbor sites 
is also reported: solid and dashed lines denote positive and 
negative pairing amplitudes, respectively. 



to ground-state properties, here we can also afford calcu- 
lations for the important low-energy excited states. We 
find that the lowest triplet excitation has different quan- 
tum numbers for J'/J > 0.775, where q = (37r/4, 0), 
and for J'/J < 0.775, where q = (tt, 7r/-\/3). Remark- 
ably, in a wide regime, the low-energy spectrum shows a 
clear ID character with two (almost) degenerate triplet 
excitations with q^ — tt, namely q — (7r,7r/V3) and 
q = (tt, 7r/3V3), see Fig. [T] 

Let us now move to a detailed study of the Hamil- 
tonian llj by using a variational wave function (WF) 
approach that includes both magnetic correlations and 
electronic pairing. In the original RVB approach the vari- 
ational WF can be obtained by applying the Gutzwiller 
projector Vg that completely suppress doubly occupied 
sites to the ground state of a mean-field BCS Hamilto- 
nian. [3| Within the same variational approach, a mag- 
netic state is obtained by adding an external field in the 
BCS Hamiltonian and considering a suitable long-range 
spin Jastrow factor that introduces the correct spin- wave 
fluctuations, i.e., J, = exp (l/2 ^-^^f . [H The 

ground state of the mean-field Hamiltonian containing 
both electronic pairing and magnetism can be written in 
terms of a generalized complex pairing function f^y ' 
that contains both singlet and triplet components, so 
that the full variational WF is given by: 



23| 



|vI')=J,7'Gexp<;- >^ f:rcl,.A,^A\0}, (2) 



where c|^. creates an electron with spin ct^ on the site 



i. At present, all variational approaches on the lattice 
have optimized the WF by considering few short-range 
parameters of the BCS Hamiltonian (e.g., the BCS pair- 
ing and/or hopping amplitudes), implying a long-range 
pairing function. Here, we generalize this variational 
approach without defining the mean-field BCS Hamil- 
tonian. Instead, we directly optimize the pairing am- 
plitude f'^j''^^ ■ This approach allows us to have more 
variational freedom and, therefore, provides a much less 
biased ansatz to the ground state. Let us now discuss 
the symmetries that we use for this quantity. First of all, 
we consider independent {ai,aj) values for (t,T)) (iji)j 
(t,i), and (|,t) amplitudes. Then, in order to take into 
account magnetic correlations, we consider two different 
possibilities. The first one has a three-sublattice symme- 
try (suitable to the 120° order) and the second one has 
antiparallel spins along ID chains and with two indepen- 
dent chains with different magnetic moment (suitable to 
describe the magnetic order for J'/ J <^1). Despite these 
limitations, the correlated WF of Eq. optimized in 
presence of the Jastrow factor J7s , may show clear incom- 
mensurate spin-spin correlations, demonstrating that our 
approach is highly fiexible and allows us to describe non- 
trivial spin correlations. In summary, for each bond and 
each spin case, we have three (two) independent com- 
plex numbers for the first (second) case. Finally, a 2 x 1 
structure that implies an extra sign factor (-1-1 or — 1) 
for each fij'^^ is considered [l3| for the case of three 
sublattices, see Fig. [2] for the sign convention of nearest- 
neighbor sites. Periodic or antiperiodic boundary condi- 
tions on fpj''^^ are chosen, depending on L. Although 
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FIG. 3: (Color on-line) Energy per site for the 12 x 12 (upper 
panel) and the 18 x 18 (bottom panel) for the WF with three- 
sublattice magnetization and 2x1 structure (red squares) 
and with two decoupled chains (blue triangles). See text for 
a detailed description of the WFs. 



the WF breaks the spin SU(2) symmetry, the actual 
value of the total spin square (5^) is as small as 0.07 
for the 18 x 18 cluster. The fundamental ingredients 
are the presence of the 2x1 structure for the signs of 
the pairing, relevant for J' ^ J, and the direct opti- 
mization of /■^■' ^ (containing both singlet and triplet 
components), which is afforded here for the first time. 
In the isotropic case, we obtain excellent results, which 
give an energy per site E/J = —0.5470(1) in the ther- 
modynamic limit, very close to our estimation of the ex- 
act value E/ J = —0.551(1) (which is extracted with the 
variance extrapolation of WFs with zero and one Lanc- 
zos step [iSl) and much lower than previous estimates 
S/J- -0.53. 



In Fig. [21 we report the accuracy of the two WFs (with 
three-sublattice or two-chain structure) for the 6x6 lat- 
tice. The full optimization of the pairing function allows 
us to reach a very good accuracy in the whole range of 
our interest and, in particular, for J'/ J < 0.5. We notice 
that the level crossing present in exact calculations (see 
Fig. [1]) is also present in the energy of the two variational 
WFs, although it is shifted to J'/J '-^ 0.6. In the case 
of a first-order transition, there is a macroscopic energy 
difference between the stable and the unstable states in 
both regions across the transition point. However, by in- 
creasing the system size, we observe that the two energies 
merge for small frustrating ratios, namely J'/J < 0.6, see 
Fig. [3l This indicates that the transition becomes con- 
tinuous in the thermodynamic limit. The tiny energy 
difference between the two WFs for J'/J < 0.6 suggests 
an effective chain decoupling. Indeed, the two variational 
WFs are compatible with a continuous transition: at the 
critical point, the two states coincide and have vanishing 
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FIG. 4: (Color non-line) Static spin-spin correlations ^(g) for 
different lattice sizes N = L x L and frustrating ratios J'/J. 
A darker color indicates a bigger S{q). 
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FIG. 5: (Color on-line) Lower panel: size scaling of the mag- 
netic order parameter for different values of J'/J. Stars indi- 
cate exact results for A'^ = 36 and lines are fits. Upper panel: 
position of the peak q = ((3,0) for J'/J = 0.8. 



inter-chain pairing amplitude. 27 1 

Let us now move to the magnetic properties that can 
be assessed by the static spin-spin correlations: 



Siq) 



(3) 



In Fig. m we show the results of S{q) for three typical 
values of J'/J and three lattice sizes. In the isotropic 
case, we found that S{q) has huge peaks at the corner 
of the Brillouin zone. The size scaling of = S{Q)/N 
with Q = (47r/3,0) indicates a three-sublattice magnetic 
order, see Fig. [5l In the thermodynamic limit, we find 
rri^ ~ 0.035, which is larger than ~ 0.02 found in 
previous works (within the present definition), [^, [l3| 
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FIG. 6: (Color on-line) Phase diagram of the anisotropic tri- 
angular lattice, as obtained by our variational approach. 



showing that our approach favors magnetic phases over 
spin liquids. Despite the fact that the magnetic moment 
is considerably overestimated, the WF captures correct 
qualitative features. For J'/ J = 0.9, we still obtain a 
finite value of m^, which is very close to the one found in 
the isotropic point. In this case, the peak of S{q) stays 
at Q = (47r/3, 0), very close to the estimation given in 
Ref. Moreover, another state can be stabilized, with 
incommensurate Q but slightly higher energy. These 
facts indicate that the true incommensurability could be 
very small and it is not detectable with the available sizes. 
On the other hand, the size scaling at J' / J = 0.8 clearly 
indicates that — > in the thermodynamic limit. Here, 
incommensurate spin correlations are found (see Fig. [4]), 
demonstrating the flexibility of the variational WF. Fur- 
thermore, for J' / J < 0.6, the spin-spin correlations dis- 
play an almost ID character: S{q) does only depend upon 
Qx, whereas it has a flat behavior as a function of qy, see 
Fig. [31 In this regime, the triplet components of the 
pairing amplitude are irrelevant, and we get a perfect 
RVB singlet state. Although we cannot exclude a tiny 
(incommensurate) magnetic order, as it was pointed out 
in Ref. 0, our calculations highlight the fact that the 
physical properties in the weakly coupled regime, i.e., 
J'/ J < 0.6, can be effectively represented as a ID spin 
liquid down to very low energies (temperatures). On the 
other hand, for J'/J > 0.6 triplet components become 
fundamental to describe magnetic fluctuations. At the 
same time, for J'/J > 0.6 also the 2x1 structure of the 
pairing turns out to be important to gain energy, indicat- 
ing a (second-order) transition between two spin liquids: 
one connected to the ID case, having all equivalent sites, 
and another one, having a 2 x 1 structure in the pairing 
function. No dimer order is found in the whole regime of 
frustration < J'/J < 1. 

In summary, by using an improved variational ap- 
proach, we have given strong evidence that a gapless spin 
liquid with negligible inter-chain coupling at low energy 
is stable over a wide region of the anisotropic triangular 
lattice. The complete phase diagram of this model (see 
Fig. [6]) can be worked out by considering both singlet and 
triplet pairing amplitudes that may give rise to magnetic 
order as well as incommensurate spin fluctuations. Our 
approach highlights the possibility to have two contin- 
uous transitions: a first one between two spin liquids 
and another one from a 2D spin liquid and a magnetic 



phase. Close to the isotropic region, triplet correlations 
are particularly important, without necessarily implying 
magnetic order. We finally remark that the existence of 
a non-magnetic state with explicit triplet pairing would 
naturally lead to a finite susceptibility at zero tempera- 
ture and the onset of triplet superconductivity, in agree- 
ment with NMR experiments in k-(ET)2Cu2(CN)3. Q 
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